LMU-07/96 
MPI/PhT/96-055 
hep-ph/9607255 
July 1996 



Dispersion Relations in Loop Calculations' 

- - - , Bernd a. Kniehl'^ 

^ : 

G^ . Institut fiir Theoretische Physik, Ludwig-Maximilians-Universitat, 



00 



TheresienstraBe 37, 80333 Munich, Germany 



Abstract 



, These lecture notes give a pedagogical introduction to the use of dispersion re- 

jy-^ I lations in loop calculations. We first derive dispersion relations which allow us to 

CN ' recover the real part of a physical amplitude from the knowledge of its absorptive 

. part along the branch cut. In perturbative calculations, the latter may be con- 

\^ [ structed by means of Cutkosky's rule, which is briefiy discussed. For illustration, 

G^ ' we apply this procedure at one loop to the photon vacuum-polarization function 

induced by leptons as well as to the 7// vertex form factor generated by the ex- 
Qh! change of a massive vector boson between the two fermion legs. We also show 

pi il how the hadronic contribution to the photon vacuum polarization may be extracted 

from the total cross section of hadron production in e~^e~ annihilation measured as 
a function of energy. Finally, we outline the application of dispersive techniques at 
the two-loop level, considering as an example the bosonic decay width of a high-mass 
^ ■ Higgs boson. 



1 Introduction 



Dispersion relations (DR's) provide a powerful tool for calculating higher-order radiative 
corrections. To evaluate the matrix element, 7}j, which describes the transition from some 
initial state, to some final state, |/), via one or more loops, one can, in principle, adopt 
the following two-step procedure. In the first step, one constructs Im7}j for arbitrary 
invariant mass, s = pf, hj means of Cutkosky's rule [Q, which is a corollary of S'-matrix 
unitarity. In the second step, appealing to analyticity, one derives Re7}j by integrating 
Im Tfi over s according to a suitable DR. 

Dispersive techniques offer both technical and physical advantages. Within perturba- 
tion theory, they allow us to reduce two-loop calculations to standard one-loop problems 
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plus phase-space and DR integrations, which can sometimes be performed analytically 
even if massive particles are involved 0, ^. This procedure can also be iterated to tackle 
three-loop problems [Q. On the other hand, dispersive methods can often be applied 
where perturbation theory is unreliable. To this end, one exploits the fact that, by virtue 
of the optical theorem, the imaginary parts of the loop amplitudes are related to total 
cross sections, which may be determined experimentally as a function of s. Perhaps, the 
best-known example of this kind in electroweak physics is the estimation of the light- 
quark contributions to the photon vacuum polarization — and thus to a(M|) — based on 
experimental data of a{e~^e~ hadrons) This type of analysis may be extended both 
to higher orders in QED and to a broader class of electroweak parameters 0. 

These lecture notes are organized as follows. In Section 2, we derive DR's appropriate 
for physical amplitudes. Cutkosky's rule is introduced in Section 3. As an elementary 
application, we calculate, in Section 4, the leptonic contribution to the photon vacuum 
polarization to lowest order in perturbation theory, and relate its hadronic contribution to 
the total cross section of e'^e~ hadrons. In Section 5, we derive, from a Ward identity, 
a subtraction prescription for general vacuum polarizations. In Section 6, we evaluate, via 
a DR, the 7// vertex form factor generated by the exchange of a massive vector boson 
between the two fermion legs. In Section 7, we outline the application of DR's at the 
two-loop level, considering the bosonic decay width of a high-mass Higgs boson. 

2 Dispersion relations 

In elementary particle physics, we often encounter form factors {i.e., functions of g^, where 
q is some transferred four-momentum) which are real-valued for below some threshold, 
M^, and exhibit a branch cut for > M^. We shall discuss various examples of form 
factors in Sections 4-7. In order to benefit from the powerful theorems available from the 
theory of complex functions, it is necessary to allow for to be complex, although this 
may contradict the naive physical intuition. 
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Figure 1: Contour C of Eq. in the complex-s plane. 
Let us then consider a complex- valued function, F{s), of complex argument, s, and 
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assume that: (1) F{s) is real for real s < M^; (2) F{s) has a branch cut for real s > M^; 
(3) F{s) is analytic for complex s (except along the branch cut). As usual, we fix the sign 
of the absorptive (imaginary) part of F along the branch cut by 

F{s + ie) =ReF{s) +ilmF{s), (1) 

where e > is infinitesimal. By Schwartz' reflection principle, we then have 

F{s + ie) - F{s-ie) =2ilmF{s). (2) 

Since F is analytic at each point within contour C depicted in Fig. |I], we may apply 
Cauchy's theorem to find 
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27rz Jc s — q 

1 ( r'-\,ni±M^iii^^j ,,ns) 



2'!Ti\jR42" s — 7|s|=A2 s — q^ ^ 

1 ImF(s) 1 / ^ F{s) 

77 Jm'^ s — q^ — le Zm j|s|=a2 s — q^ 

where we have employed Eq. (|]) in the last step. Suppose that we only know ImF along 
the branch cut and wish to evaluate F at some point q^. Then, Eq. (^) is not useful for 
our purposes, since F also appears on the right-hand side, under the integral along the 
circle. Thus, our aim is to somehow get rid of the latter integral. If 

lim / ds = 0, (4) 

A2^oo/|s|=A2 S- q^ 

then we obtain the unsubtracted DR 

^/ 9n ^ f°° , ImF(s) 

F{q') = - ds (5) 

TT JM2 s — q"^ — le 

This means that F can be reconstructed at any point q^ from the knowledge of its ab- 
sorptive part along the branch cut. In particular, the dispersive (real) part of F may be 
evaluated from 



ReF{q') = -V ds ^, (6) 

TT JaP s — q^ 

where V denotes the principal value. 

Equation (^ is not in general satisfied. It is then useful to subtract from Eq. (^) its 
value at some real point q^ < M^, 

2, a2_^2 A2 j^^^^) ,q^-ql[ , F{s) 



Fiq') = Fiqo) + - — - V ■ + V-^ f ds — ^. (7) 

TT JkP s - q^ s - q^ - le 2m J\s\=K^ {s - q^){s - q^) 
If the last term in Eq. (0) vanishes for — oo, then we have the once-subtracted DR 

- ds ImF(s) 



M'^ s — ql s — q^ — ie 

Otherwise, further subtractions will be necessary. For the use of DR's in connection with 
dimensional regularization, we refer to Ref. 
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3 Cutkosky's rule 



In the previous section, we explained how to obtain the dispersive part of a form factor 
from its absorptive part. Here, we outhne a convenient method how to evaluate the 
absorptive part within perturbation theory. 

Decomposing the scattering matrix as 5 = 1 + iT, where T is the transition matrix, 
we obtain 

- z(T - T^) = T^T (9) 

from the unitarity property S'^S = 1. Since four- momentum is conserved in the transition 
from some initial state \i) to some final state |/), we may always write 

(/|TK) = (27rr5W(P;-P,)r;,. (10) 

Consequently, 

(/|rt|.) = {^\T\fr = (27r)V^)(P^ - P,)T*. (11) 
Inserting a complete set of intermediate states \n), we find 

(/|TtT|z) = ^(/|Tt|n)(n|TK) 

n 

= (27r)V^)(P^ - P,) 5:(27r)V^)(P„ - P.)r„*^r„,. (12) 

n 

Using Eqs. ([T0|)-(p!2D in connection with Eq. (|^) and peeling off the overall delta function, 
we obtain Cutkosky's rule, 

- ^(^/. - Vf) = E(2vr)V^^(Pn - P^)r:f%^^, (13) 

n 

where it is understood that the sum runs over all kinematically allowed intermediate states 
and includes phase-space integrations and spin summations. Appealing to time-reversal 
invariance, we may identify the left-hand side of Eq. ([T3| ) with 2Im7}j. In summary, we 
may construct the absorptive part of a loop diagram according to the following recipe: 
(1) cut the loop diagram in all kinematically possible ways into two pieces so that one 
of them is connected to \i) and the other one to |/), where cut lines correspond to real 
particles; (2) stitch each pair of pieces together by summing over the spins of the real 
particles and integrating over their phase space; (3) sum over all cuts. 

4 Photon vacuum polarization 

One of the most straightforward applications of Cutkosky's rule and DR's is to evaluate 
the one-loop photon vacuum polarization induced by fermions. To avoid possible compli- 
cations due to large nonperturbative QCD corrections, we start by considering leptons. 
For the sake of generality, we keep the electric charge, Q, and the number of colours, Nc, 
arbitrary. 
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Figure 2: Application of Cutkosky's rule to the lepton- induced photon vacuum 
polarization. 

4.1 Leptonic contribution 

Given the QED interaction Lagrangian, £/ = —eQipj^ip, we wish to compute Tfi depicted 
in Fig. ^. We start from the cut amplitudes, 

i%i = u{k){-ieQ)-i^v{k), 

i%f = u{k){-ieQ)Yv{k). (14) 

Summation over spins yields 



spins 



AN^e^Q^ ( fc^r + - ^g^"" ) , (15) 



where m is the fermion mass and s = = 2{k ■ k + m^). Defining dk = {d^k/ (2TT)^2k' 
Eq. (13) takes the from 

21mTfi = J dkdk {2nY5^^\p -k-k)J2 Vf^m 

spins 



= 2N,aQ^sll-^T^\ (16) 
where a = (e^/47r) is Sommerfeld's fine-structure constant and 

T^^^ = y ^ (^k'^k'' + k^'k" - ^g'"'^ . (17) 

Upon integration, p is the only four-momentum left, so that we can make the ansatz 
rpnu _ jipfj-pi^ _|_ Bg'^'^. In order to determine A and B, we form the Lorentz scalars 



p^p^Tf'" = sisA + B) = 
g^^T'"' = sA + AB = -s -2m^. 
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We so find B 



-A/s 



-s — 2m^, so that IiiiTa- 



Nr. 



Im7r(s) = — aQ 




p'^p'')lm.7[{s), where 



(19) 



is the absorptive part of the photon vacuum polarization. Notice that the dot product of 
Tfi with p vanishes in comphance with electromagnetic gauge invariance. Using the once- 
subtracted DR d^) with gQ = 0, we obtain the renormalized photon vacuum polarization 



as 



vr s 



TT[S) - 
S rco 



7r(0) 

ds' Imvrfs'l 



TT J4m2 S S 
Nr. „2 " 



s — te 



(20) 



where 



fir) = -{2 



1 



arsinh 



- 5 1 

-r + - + -: 
6 r 



(21) 



r 

r/ V r 

appropriate for r < 0. This agrees with the well-known result found in dimensional 
regularization Representations of / appropriate for < r < 1 and r > 1 emerge from 
Eq. ( pll) through analytic continuation. Specifically, we have 




We verify that, for r > 1, Eq. ( [TP| ) is recovered from Eq. (^0^. The expansions of / for 
Irl -C 1 and r ^ 1 read 



(23) 



f{r) = -r + 0{r'), 
5 

Re f(r) = - ln(4r) + - + O (- 

3 Vr 



respectively. From Eq. (|23D , we conclude that heavy fermions, with mass m ^ y|s|/2, 

decouple from QED [^, while light fermions, with m <^ -y/jsj/2, generate large logarith- 
mic corrections. The latter point creates a principal problem for the estimation of the 
hadronic contribution to vr. The evaluation of Eq. ( PD| ) using the poorly known light-quark 
masses, rrig, would suffer from large uncertainties proportional to drriq/mq. In addition, 
there would be large nonperturbative QCD corrections in connection with the subtraction 
term 7r(0). In the next section, we discuss an elegant way to circumvent this problem. 
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4.2 Hadronic contribution 



Let us consider the creation of a quark pair by e^e annihilation via a virtual photon, 
e~(Z)e+(Z) — > 7*(p) q{k)q{k). The corresponding T-matrix element reads 

tT = v{l){-ieQ,)^^u{l)^^i{Tnir, (24) 

s + le 

where s = and i{TniY is given in Eq. ([l^). Taking into account the spin average 
and the flux factor, we evaluate the total cross section as 



a{s) = \y j ^'^^'^ {2Tif5^^\l + l-k-k)Y. 1^1' 

^ spins 

tr /7. j d'~kd'h {27r)^6^'\p -k-k)Y. (T„*)^(r„ 



spins 

mj 

spins 

""'^^ 4 + l^U - '-g^u) {-2){sg^^ -p^^p-) Im-nis) 



8s3 



8s3 



Im7r(s), (25) 



s 

where we have exploited Eq. ([16|) . Substituting Eq. (^Sj) into Eq. (^), we obtain 

7r(s) = r rfs'^^^. (26) 

Equation (^) allows us to estimate the hadronic contribution to tt from the total cross 
section of hadron production by e^e~ annihilation measured as a function of s. A recent 
analysis § has yielded -7^{M'^z)\hadrons = 0.0280 ± 0.0007. 



5 Subtraction prescription for general vacuum polar- 
izations 

Vacuum polarizations have mass dimension two, so that the unsubtracted DR (^ in 
general leads to ultraviolet divergences quadratic in the cutoff A, which violate the Ward 
identities of the theory and are not removed by renormalization. It is therefore necessary 
to use a subtraction. As we have seen in Section 3, in the case of QED, the naive 
subtraction of Eq. (|^) with = leads to the correct physical result. In the presence 
of unconserved currents, as in the Standard Model, the situation is more complicated Q . 



In the following, we discuss a suitable subtraction prescription |]TT . 

Starting from the general interaction Lagrangian Cj = —B^J^ + h.c, where is 
some vector field and is the associated current, it is straightforward to derive the 
vacuum-polarization tensor, which, by convention, differs from the T-matrix element of 
B'^{p) — > B^Ip) by a minus sign, as 

W{p) = -i j d^xe'P-''{Q\TJ^'{x).r\Q)\Q), (27) 
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where T denotes the time-ordered product. By Lorentz covariance, li^'^ has the decom- 
position 

n^'^(p) = n(s)^^'^ + a(s)pV, (28) 

where s = p^. Integrating by parts, we obtain from Eq. ( P?] ) 



p^W{j)) = j (i^xe*f-^(0|r9^J^(x)J'^t(0)|0) 
= p^A(s), 



(29) 



where the last step follows from Lorentz covariance. On the other hand, from Eq. (p8|) 
we get p^Ii^'^ip) = p'^[Il{s) + sX{s)]. Consequently, we have 



n(s) = A{s)-sX{s) 



(30) 



By its definition (^), 11 has mass dimension two, so that its evaluation from Imll 
via the unsubtracted DR (^ would be quadratically divergent. However, Eq. (|30|) relates 
n to quantities for which unsubtracted DR's are only logarithmically divergent. This 
is obvious for A, which is dimensionless. In the case of A, this may be understood by 
observing that is softly broken by mass terms and that one power of the external 
four- momentum is extracted in Eq. (p9l). Writing unsubtracted DR's for A and A, we find 



His) = - 



A2 



ds' 



TT JM"^ S 
1 rA-^ 



ds 



TT JM2 



■ [ImA(s') - slmA(s')] 

S Ij^ 

Imn(s') 



s' — s — le 



+ ImA(s') 



(31) 



where is the lowest threshold. Detailed inspection at 0{a.) and 0{aas) [|r^ reveals 
that the logarithmic divergences of Eq. (|3TD exhibit a similar structure as the poles in 
e = 2 — n/2 found in n-dimensional regularization, so that both methods lead to the same 
physical results. 



6 Vertex correction 

Cutkosky's rule (p!3D is also applicable if the initial and final states are different. As 
an example, we now consider the interaction between a massless-fermion current and a 
neutral vector boson, B, with mass M characterized by the Lagrangian Cj = —gip^ip, 
where g is the coupling constant, and evaluate the Bff vertex form factor induced by the 
exchange of B between the two fermion legs as depicted in Fig. |^. 
The cut amplitudes read 

iTni = u{k - q){-ig)'-i^v{k + g), 

i%S = -u{k - q){-tg)Yu{k) vCk)i-^9hMk + q) ^ . . (32) 
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Figure 3: Application of Cutkosky's rule to the Bff vertex form factor 
induced by the exchange of B between the two fermion legs. 

Notice the extra minus sign of 7^/. The spin summation yields 



9' 



spins 



u{k)T''v{k), 



(33) 



where 



r 



;— 4)'~i^{^+ ifj'ju- Anticipating the multiplication with the Bff tree-level 



amplitude, we can bring into the form = N'j^, where 



iV = --tr ^T^h, 
g 

= -k ■ {k + q)k ■ {k — q) 



(34) 



and s = p^. The integration over the phase space of the cut particles, with four-momenta 
k^ = k + q and k_ = k — q, is most conveniently carried out in the centre-of-mass frame 
using the variables p = k^ + k^ and r = k^ — k_. Then, we have = — s(l + z)/2 and 
N = s(l — z)'^/2, where z is the cosine of the angle between r and k. Cutkosky's rule (p!3| ) 
now takes the form 

ImT^, = i 1 d~k+d~k. {2nY5^'\p - fc+ - fc„) ^ r„=)T„, 



where 



spins 



-gu{k)^^v{k)^\mF{s 



(35) 



ImF(s) 



2 

TT 



-1 l + z + 2/x 
2^1 + ^ ln(l + a;) -3 



(36) 
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and X = s/M'^. Finally, we obtain the renormalized vertex function through the once- 
subtracted DR (H) as 

=F(s) -F(0) 

s ds' ImFis') 



IT Jo s' s' — s ~ te 

2 



= 2 (l + ^) [Li2(l + x)- C(2)] + (3 + ^) ln(-a:) - ^ - | (37) 

where Li2(-2) = — Jq dt \n(l — zt)/t is the dilogarithm. This agrees with the corresponding 
result found in dimensional regularization |]T2[. Notice that, in Eq. (P7D, x = {s + ie)/M'^ 
comes with an infinitesimal imaginary part. The dispersive and absorptive parts of F for 
X > are conveniently separated by using the relations 

ln(— a;) = Ina; — in, 
Li2(l + x)- C(2) = - Li2(-x) - In(-x) ln(l + x). (38) 

We so recover Eq. (PB|), which serves as a useful check. 



7 Two-loop application 




Figure 4: Massive scalar two-loop three-point diagram, which contributes to 
the Higgs-boson decay to a pair of intermediate bosons, a, b, and c denote the 
squared masses of the respective loop particles. 

In this final section, we illustrate the usefulness of DR's beyond one loop, considering 
the massive scalar two-loop three-point diagram depicted in Fig. ^ The idea is to write 
the one-loop two-point subdiagram contained in that diagram as a once-subtracted DR 
and to interchange the loop and DR integrations, so as to reduce the two-loop problem 
at hand to a one-loop one with a subsequent DR integration |I3| . 
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In dimensional regularization, the scalar one-loop two-point integral is given by 

1 



B{s, a, h) 



477 / 



(47r)2 

i 

(4^ 



(27r)" (g2 _ a ie) [{q + - b + ie] 
1 dx 

1 



+ fis) + Oie] 



(39) 



where 7 is the Euler-Mascheroni constant, F is the gamma function, n = 4 — 2e is the 
dimensionality of space time, /i is the 't Hooft mass scale introduced to keep the coupling 
constants dimensionless, s = p"^, X = [(1 — x)a + xb — x{l — x)s — ie] //x^, and /(s) = 
— Jq dx \nX. The peculiar form of the prefactor in Eq. ( pQ]) is to suppress the appearance 
of the familiar terms proportional to 7 — ln(47r) in the expressions. The absorptive part 
of / is 



Im/(s) = dx9(^x{l — x)s — (1 — x)a — xb 



e(s-(v^ + ^/b 



(40) 



where A(s, a, b) = s"^ + + b"^ — 2{sa + ab + bs) is the Kallen function. Using Eq. (^, we 
then find 



B{s, a, b) — B{c, a, b) 



(477) 

i 

(4^ 
+ Oie 



-e^^Til + e) [/(.)- /(c) + 0(5)] 



1 

71 Jo 



da Im/(cr) 



a — s — le a — c — le 



(41) 



Consequently, the subdiagram in Fig. | consisting of the bubble and the adjacent propa- 
gator with four-momentum q and squared mass c may be written as |T^ 



B{q'^,a,b) B{c,a,b) 



q'^ — c + ie q'^ — c + ie (Air)'^ 
+ 0{e). 



e^^F 1 + e 



da 



\/A(cr,a,6) 



(Va+Vb) a — c — it 



a 



a + ie 



(42) 



If Eq. ( ^2]) is inserted in the expression for the one-loop seed diagram in Fig. ^ the first 
term turns into a product of two one-loop diagrams, which contains all divergences. In 
the second term, we may interchange the DR and loop integrations and are left with a 
finite dispersion integral, which may be solved analytically. 
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